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O Abstract 

We study Dirac spinors in Bianchi type-I cosmological models, within the 
framework of torsional /(i?)-gravity. We find four types of results: the re- 
sulting dynamic behavior of the universe depends on the particular choice of 
function f(R)] some f(R) models do not isotropize and have no Einstein limit, 
so that they have no physical significance, whereas for other f(R) models 
^ H isotropization and Einsteinization occur, and so they arc physically accept- 

5— i able, suggesting that phenomenological arguments may select f(R) models 

1 ^9 that are physically meaningful; the singularity problem can be avoided, due 

to the presence of torsion; the general conservation laws holding for f(R)- 
gravity with torsion ensure the preservation of the Hamiltonian constraint, so 
<*" proving that the initial value problem is well-formulated for these models. 

r^j- Keywords: Dirac spinors, Bianchi-I models, /(i?)-cosmology, torsion 

O PACS number: 04.50.Kd, 98.80.Jk, 04.62. +v 

2 1 Introduction 

The need to enlarge or revise General Relativity (GR) arises from some open prob- 
lems in modern physics at both ultraviolet and infrared scales. For instance, at 
astrophysical and cosmological scales the predictions of GR agree with observa- 
tions only after the theoretical assumption regarding the existence of the so-called 
dark matter and dark energy. But at fundamental level, up to now there are no 
experimental evidences that such unknown forms of matter and energy really exist. 
This fact together with other shortcomings of Einstein's theory represents signals 
of a possible breakdown in our understanding of gravity, and the opportunity of 
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developing extended or alternative theories of gravity is to be seriously considered. 
Actually, several extensions of GR have been formulated in the last thirty years; 
among these, /(-R)-gravity is one of the most direct and simplest: it consists in 
relaxing the hypothesis that the gravitational Lagrangian be a linear function of 
the Ricci scalar R, thus allowing a dependence on R more general than that of the 
Einstein-Hilbert action. In the last years, /(-R)-gravity has received great atten- 
tion due to the fact that it seems to account for both cosmic speed-up and missing 
matter at cosmological and astrophysical scales, respectively [H [21 [3] . 
At the same time, another way to extend GR consists in considering torsion, includ- 
ing it among the geometrical properties of space-time. Torsion was first introduced 
in the geometrical background by Cartan; then Sciama and Kibble embodied it 
into the framework of Einstein gravity following the idea that energy is the source 
of curvature in the same way in which spin is the source of torsion. The resulting 
Einstein-Cartan-Sciama-Kibble (ECSK) theory has been the first generalization of 
GR trying to take the spin of elementary fields into account and it is still one of 
the most serious attempts in this direction [U [5] . 

According to this line of thinking, /(fi)-gravity with torsion represents one of 
the simplest extensions of the ECSK theory, just as purely metric /(ii)-gravity is 
with respect to GR. Indeed, the basic idea is again that of replacing the Einstein- 
Hilbert Lagrangian with a non-linear function. An important consequence of the 
non-linearity of the gravitational Lagrangian is that we have non-vanishing torsion 
even without spin, as long as the stress-energy trace is not constant [6], [U |9] . It 
is known that torsion may give rise to singularity-free and accelerated cosmologi- 
cal models [10]; just arising from the non-linearity of the gravitational Lagrangian 
function f(R), a propagating torsion could amplify this effect and this point cer- 
tainly deserves further attention. 

Besides, in a recent paper [llj we have studied in detail /(i?)-gravity with torsion 
in presence of Dirac fields, in order to better exploit the coupling between torsion 
and spin within the framework of /(i?)-theories. 

In the present work we explore possible cosmological applications of the theory 
proposed in More precisely, we study Dirac fields in Bianchi type-I (BI) f(R)- 
cosmological models. BI is the simplest anisotropic space-time that generalizes the 
spatially flat Friedmann-Lemai'tre-Robertson- Walker (FLRW) universe. The dif- 
ference with respect to the FLRW flat model is that in a BI universe the scale 
factors of spatial directions can differ from each other, thus yielding anisotropy. 
Though FLRW models provide very accurate descriptions of the universe observed 
today, the latter could have undergone an anisotropic phase in its early epoch; 
for instance, the detected anisotropies in the cosmic microwave background could 
represent evidences in this direction. On the other hand, non trivial Dirac fields 
are seen to generate spin-torsion which is incompatible with spatial isotropy |12j : 
when dealing with Dirac fields in presence of torsion, anisotropic models are then 
a forced choice. Of course, for consistency with the observed data at present day, 
the initial anisotropic phase have to undergo a subsequent isotropization process 
as the Dirac fields energy and spin decrease and universe expands. This dynamical 
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behavior can represent a selection rule for viable functions f(R). 
The layout of the paper is the following. In Sections 2, we review some gener- 
alities about /(.R)-gravity with torsion in presence of Dirac fields. In Section 3, 
we introduce BI cosmological models in torsional /(-R)-gravity still in presence of 
Dirac fields and, for completeness, perfect fluid. Following the line traced in pre- 
vious papers [131 Q2]> we derive a solution of the generalized Dirac equations, 
holding independently of the assumed model f(R)', we integrate the Einstein-like 
equations in the simplest case f(R) = R and discuss an approximate solution for 
f(R) = R + aR 2 . Section 4 is devoted to discuss the field equations in the Einstein 
frame; we show that, through a conformal transformation from the Jordan to the 
Einstein frame, we can put the Einstein-like equations in a quadratically integrable 
form; this result turns out to be useful for a qualitative analysis of the solutions 
or numerical integration even when explicit calculations for analytical integration 
are difficult; two examples are given. Finally, in Section 5, we face a crucial as- 
pect in dealing with Einstein-like equations, namely the preservation in time of 
the Hamiltonian constraint (momentum constraints are here trivial identities); in 
GR this result is generally ensured by the Bianchi identities and the matter field 
equations, but here it is not a priori expected and it has to be proved; anyway, by 
working out the general conservation laws proved in |llj . we show that classical 
results by Bruhat [IB] apply also in the present case, thus ensuring the required 
time-preservation of the Hamiltonian constraint. 



2 Dirac fields in /(.R)-gravity with torsion 

In previous papers [61 [T7], /(-R)-gravity with torsion has been formulated 
within both the metric-affine and the tetrad-affine framework. According to the 
basic paradigm of /(i?) -gravity, the gravitational Lagrangian of the theory is as- 
sumed to be a real function f(R), where R is here the Ricci curvature scalar 
written in terms of a metric g and a g-compatible connection T, or equivalently a 
tetrad field e and a spin-connection oj. The pairs (g, T) and (e,w) represent the 
gravitational dynamical fields of the theory in the metric-affine and tetrad-affine 
approaches respectively. The field equations turn out to be 
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f'(R)R ij - ^f(R)g l3 = X, 



(2.1a) 



f'{R) 



df'(R) 
dxP 



' pq 



+ 5, 



(2.1b) 



where Rij and T { - h are the Ricci and Torsion tensors, while and S^ h are the 



stress-energy and spin density tensors of the matter fields. From (2.1b), it is seen 



that we can have non-vanishing torsion even in absence of spin density. 

Making use of the Bianchi identities, it is possible to derive the conservation laws 
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(2.2a) 
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V h S ijh + T h S ljh + Z ij - TP* = (2.2b) 

under which the stress-energy and spin density tensors of the matter fields must 
undergo once the matter field equations are assigned [11]. In equations (2.2) the 
symbols Vj denote covariant derivative with respect to the dynamical connection 
r while R st%:) is the curvature tensor of V. Indices are lowered and raised by the 
metric gij. 

In this paper we deal with /(i?)-gravity coupled with Dirac fields and possibly a 
perfect fluid (without spin). Denoting by 7 M (fi = 0,1,2,3) Dirac matrices, we 
introduce the notation P = where ef (efel = by and e^ej = <5j) indicate a 
tetrad field associated with a metric g^. Moreover, setting S^ u := | [7^,71/], we 
denote the covariant derivative of the Dirac field ip by D^ip = ^ + uj, i ^ v S^ u ip and 
Diip = — ipLOi ^ S^ u , where u i ^ u is a spin connection. Equivalently, we can put 
Diip = & ~ V-iip and = S- + ip^ where 



^■■=-\ 9 j h (rJ-e^ k )r h r k (2.3) 

T ik being the coefficients of a linear connection V. The relation between linear 
and spin connection is given by 

T ii h = Ui \^ + ^ (2.4) 
The stress-energy tensors of the matter fields are then described as 

Eg := - (rpTiDjip - D^T^) (2.5a) 
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and 



i j 

while the spin density tensor is expressed as 
i - r„, „ -i . l 



Xij-iP+pWiUj-m, (2.5b) 



S H h = \$ { r "' *%} ^ = "4»T^Ar (^757 T ^) e£<^ (2.5c) 

In equations (2.5), p, p and Ui denote respectively the matter-energy density, 
the pressure and the four velocity of the fluid, Sij := ^[Ti,Tj] and rf a is the 
Minkowskian metric with signature (1,-1,— 1,-1). 

Following previous works [H El EJ EJ [H] , we suppose that the trace of equations 



(2.1a 



f(R)R-2f(R)=X (2.6) 

gives rise to an invertible relation between the Ricci scalar curvature R and the 
trace £ of the stress-energy tensor. Also, we assume that f(R) ^ kR 2 (the case 
f(R) = kR 2 is only compatible with the condition E = 0). Under the assumed 
conditions, from equation ( |2.6| ) it is possible to derive the expression of R as 
function of S, namely R = F(E). After that, introducing the scalar field 

<P ■■= f'(F(Z)) (2.7) 
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and the effective potential 



V&) : = - [^((/TV)) + tfUTKv)] 
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as well as separating the Levi-Civita contributions from the torsional ones, we can 
express the field equations (2.1a) in the Einstein-like form 

3 dp dp 



R; 



-RQij 



J ^3 J p A 



2 dx l 8x3 



^_ dtp 3 dp dp hk 



~^ h d% 9 - 



V(p)g lJ \+V h S ji h + S hi p S : 



JV 



1 

2"hq 



Si, p S q p gij 



(2.9) 



where Rij, R and V, denote respectively the Ricci tensor, the Ricci scalar curvature 
and the covariant derivative of the Levi-Civita connection, and := —^S i j h . 
In addition to this, the generalized Dirac equations for the spinor field are 



iY h D h ip + -T h T h ^ -mip = 



(2.10) 



where := T^P is the torsion vector. As it has been proved in [TT], equations 
(2.10) imply automatically the validity of the conservation laws 

1 



v,->: / "- / + l)y. n 'i 



S sU R stij = 



(2.11a) 



VhS ijh + ThS ijh + ^Dij _ ^Dji = q 



(2.11b) 



Moreover, still in |11| . it has been shown that equations (2.11b) are equivalent to 
the antisymmetric part of the Einstein-like equations (2.9). From this, we get two 
results: on one hand, we have that the significant part of the Einstein-like field 
equations (2.9) is the symmetric one; on the other hand, comparing equations 
(2.11a) with equations (2.2a), we derive the conservation law for the perfect fluid. 
Indeed, introducing the contorsion tensor 



Ki 



h 



rp h I rp h 

ij 3 i 



-■ft 



and expressing the dynamical connection as the sum 



(2.12) 



h 



ij 



13 



K 



h 



'J 



(2.13) 



where are the coefficients of the Levi-Civita connection, a direct comparison 
between equations (2.11a) with equations (2.2a) yields the identity 



+Th ^Fhj _ T j ih Z Fih = V^ Fij = 



(2.14) 



5 



which is the same conservation laws holding in general relativity. 
To conclude this preliminary section, we notice that the symmetrized part of the 
Einstein-like equations (2.9) as well as the Dirac equations (2.10) can be worked 
out as in [11] assuming the final form 



R, 



-Rg 



dp 



3 dtp dip ~ dp 3 dtp dp 



2 dx' 1 dxi 



dx 1 4 dx h dx k 



;9 hk gij 



Up 2 



and 



where 
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iT h D h ip — [(i>ip) + ^(#75^)75] tp-mip 



(2.15) 
(2.16) 



(2.17) 



Di denoting covariant derivative with respect to the Levi-Civita connection. 

3 Bianchi-I cosmological models 

Let us consider a Bianchi type I metric of the form 

ds 2 = dt 2 - a 2 (t) dx 2 - b 2 (t) dy 2 - c 2 {t) dz 2 
The tetrad field associated with the metric (3.1) is given by 



e$ = 8$, e» = a(t)5l e£ = b(t)6$, e» 
The inverse relations of (|3.3l) are expressed as 



„0 



1 



a(t) 



1 

W) 



\i = 0,1,2,3 



(3.1) 
(3.2) 



^ = 0,1,2,3 (3.3) 



The non-trivial Christoffel symbols associated to the metric (3.1) are 



f 1 
1 10 



a 
a 
ad, 



20 
f 

11 — UU, 1 22 

In this case, the matrices T' 1 = e l 7^ assume the explicit form 



_ b 
= V 
66, 



f 3 

1 30 



f 
1 33 



C 
C 
CC 



(3.4) 



r 1 



1 



r7 



1 



:7 



r7 



(3.5) 



a{t) 1 ' b(t) 1 ' c(t) 

The spinorial covariant derivative induced by the Levi-Civita connection is 

Ditp = diip - Qiijj, Di-ip = dii) + 4>tti (3.6) 
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where the spinor connection coefficients are given by 



(3.7) 



Taking equations (3.6) and (3.7) into account, it is easily seen that the Dirac 
equations (2.16) assume the form 

^ + + WV - T7T + i^TVhV] i> = ( 3 - 8a ) 

2r lb(f 

t + - imipj + [(#)7° + i(#7 5 V07 5 7°] = (3.8b) 

It Vop> 

where, borrowing from |13l I15j . we have defined r := abc. Analogously, evaluating 
the Einstein-like equations (2.15) for the metric (3.1), we get 

3 



ab be ac p 1 - 

-- + -- H = - + —mipip 

ab be ac p zip 



Up- 



+- 



1 

7 2 



r (^757^) (^757^) + 

ft _ _ V(<p) 

4 T 



(3.9a) 



6 c 6 c 
o c be 



a c ac 
- + - + — 
a c ac 



a b ab 
a b ab 



V 1 

— + -2 



p 1 



W: + 7V 2 ~ v ( v 3 + L( P ) ~ 



+ 

ip ip* 



a 4' 



6 3 



+ 



64v?- 



+ 



:(V ; 757 T V')(V , 757tV ; ) 



b 4 



c 3 



+ 



+ 



6499- 



-(^757 r V , )(V ; 757r^) 



c 4' 



+ 7^ - V 3 V 3 + -<P - V(y>) 



+ 



+ 



64^- 



:(V'757 T V ; )(V'757rV ; ) 



Sf 2 = 

sf 3 = o 

Sf 3 = 



(3.9b) 



(3.9c) 



(3.9d) 



together with the conditions 

» a6-6d = U tj)j 5 j 3 ip = (3.10a) 

c 6-6c = U V¥V^ = (3.10b) 

» ac-cd = U V5 7 5 7 2 ^ = (3.10c) 

The equations = (A = 1, 2, 3) are automatically satisfied identities. Finally, 
the conservation law (2.14) together with an equation of state of the kind p = Xp 
(A G [0, 1[) yield the last equation 



p + -(l + A)p = 
r 



(3.11) 
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which completes the whole set of field equations. The general solution of (3.11) is 
given by 

po = constant (3-12) 



P = POT 



-(1+A) 



Conditions (3.10) are constraints imposed on the metric or on the Dirac field. 



We see that there are in general three ways of satisfying these conditions: one is 
to impose constraints of purely geometrical origin by insisting that ab — ba = 0, 
ac — ca = 0, cb — be = giving an isotropic universe filled with fermionic matter 
fields, which is problematic due to the fact that it is known that Dirac fields do not 
undergo the cosmological principle [T2] ; another is to impose constraints of purely 
material origin by insisting that tfj^^ip = 0, Tpj 5, y 2 i/j = 0, ip"f 5, y 3 ijj = giving an 
anisotropic universe without fermionic torsional interactions, which we regard as 
unsatisfactory because if Dirac fields are absent then it is not clear what may then 
justify anisotropies; the last situation would be of both geometrical and material 
origin by insisting that for instance ab — ba = with ?/>7 7 -i/; = 0, tjj^j 2 ^ = 
giving a partial isotropy for only two axes with the corresponding two components 
of the spin vector vanishing, describing a universe shaped as an ellipsoid of rotation 
about the only axis along which the spin vector does not vanish. Notice that by 
insisting on the proportionality between two couples of axes we inevitably get the 
total isotropy of the 3-dimensional space. Therefore, the situation in which we 
have a = b with tfrf r f-ij) = ■ip^ 5 'y 2 'ip = is the only one that we believe to be 
entirely satisfactory, and from now on we shall work in this situation. 



Here, the Dirac and Einstein-like equations (3.8) and (3.9) can be worked out as 
in |13} I15j : for instance, through suitable combinations of (3.9) we obtain the 
equations 

d_ 
eft' 
3 



rft/ftyr) = 



d_ 
dt 



(#7 5 '0)(V ,t 7 5 V'T") = 



(3.13a) 
(3.13b) 



dt 



2m + 



Tt 



(-0 t 7 5 V'T) - 2m(iipj 5 ipT) = 



From equations (3.13) it is easy to deduce that 



(^ 7 Vvr = (W>) + (#rV) + (V>7 7» 



r 2 



K 2 



(3.13c) 
(3.13d) 

(3.14a) 
(3.14b) 



with C and K constants. We notice that in this special case the theory has an 
additional discrete symmetry given by the fact that under the discrete transforma- 
tion ip — > 7 5 7°7 1, i/; all field equations are invariant; this implies that in the Dirac 
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equation the total number of 4 components is in this case reduced to 2 components 
alone; however 2 components with complex values are equivalent to 4 components 



of real values: so the 4 equations for real fields given in (3.13) are the system of 



field equations we will have to solve. The compatibility with all constraints allows 
only three classes of spinors, each of which has a general member written in the 
following form 



4 



with constraints tan C,i tan £2 
integer, and also 



1> 



and 



!_ 

2t 



/ y /K -C cosCie iei \ 
VK + C cos C 2 e i,?1 
v^-CsinCie^ 2 
V v / ^ + CsinC 2 e i92 / 



( — l) n+1 and 8\ + 82 — $1 — $2 = for any n 
( VK -CcosCie i01 \ 



2r 








(3.15) 



\ v / ^ + CsinC 2 e i92 / 



i_ 

2r 



/ 

VKTC cos de^ 1 
v / ^-CsinC 2 e^ 2 



V 







(3.16) 



where all angular functions Ci, C 2 an d ^1, O2, ^ 2 have only temporal dependence, 
that has to be determined by plugging the spinor back into the Dirac equations. 



For the gravitational field, borrowing again from [15] , we subtract equation (3.9b) 



from equation (3.9d) obtaining 



tpr 



d_ 
dt 



+ ipr 



+ fir 



from which we derive 



(3.17) 



(3.18) 



D and X being suitable constants, which gives the evolution equation for the 
ratio of the two axes that are not constantly proportional, and therefore it can be 
considered as the evolution equation for the shape of the universe; also, multiplying 



3.9a) by 3 and adding the result to the summation of equations (3.9b), (3.9c) and 



3.9d), we get the final equation for r 



r ld 
2- + 3- 
r if 



p 3 
!- + — 

ip 2tp 



rrnpip — 5 



T If 

T if C/2 Z 



>V(<p) 



(3.19) 
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which can be considered as the evolution equation of the volume of the universe. 



Here, it is worth noticing that equation (3.9a) plays the role of a constraint on 



the initial data: thus for consistency we have to check that, if satisfied initially, 
this constraint is preserved in time. To see this point, we first observe that the 



Einstein- like equations (2.15), and thus also (3.9), can be written in the equivalent 
form 

1. 



T 



-Tgi 



(3.20) 



where 



-%3 •- 



ip J p J p z 



3 dip dp ^ dip | 3dip_dp_ hk 

+ (pV >»~ i + 4dx>>dx k9 9lj 



2 dx % dxi r 3 dx 



(3.21) 



Up- 



denotes the effective stress-energy tensor appearing on the right hand side of equa- 



tions (2.15), while T is its trace. It is then a straightforward matter to verify that 



equations (3.17) and (3.19) can be equivalently obtained by suitably combining the 



space-space equations of the set (3.20); more in detail, equation (3.17) is obtained 



by subtracting each to other the two distinct space-space equations of (3.20), while 



equation (3.19) is obtained adding together all the space-space equations of (3.20). 



As a consequence, we have that solving equations (3.17) and (3.19) amounts to 



solve all the space-space equations of the set (3.20). In addition to this, as it is 



proved in Section 5, the conservation laws (2.2) automatically imply the vanishing 



of the four-divergence of the Einstein-like equations (2.15). The two just mentioned 



facts allow to apply to the present case a result by Bruhat (see [16], Theorem 4.1, 



pag. 150) which ensures that the constraint (3.9a) is actually satisfied for all time 



Examples: f(R) = R and f(R) = R + 5aR 2 models. In the following, we will 
consider the situation given by f(R) = R + 5aR 2 as a correction in da with respect 
to the f(R) = R case, and we are going to take solutions for the Einstein-like field 
equations that are corrections in 5a of the exact solutions of the Einstein field 
equations: thus the first step is to find the exact solutions of the f(R) = R case. 
To this extent, consider for simplicity the simplest case without fluid. Then, we 



choose for the Dirac field the solution in the form given as in (3.16) by 



2r 







y/K + Cj 



-vat- 



VK - Ce 



V 



o 



mt- 



3U f d 
16 J 7 

3C r . 
16 J 



J 



(3.22) 



with constraints xptp = —, ifiip = — and ^^°ip = 0, i^j 5 ifj = 0, and with 
the solution defined up to a global unitary phase which is then irrelevant: this 
solution albeit simple is nevertheless not trivial, and it will actually lead to a 
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correspondingly simple form of the gravitational evolution equations for the shape 
and volume of the universe as 



(3.23) 



and also 



3mC 







(3.24) 



in terms of the constants m, C, X, D alone. As it is clear, it is possible to integrate 



the equation (3.24) to give the evolution of the volume of the universe as 

r = -—(b + 2(3t + t 2 ) 



(3.25) 



where b is the integration constant that contains the information about the initial 
volume of the universe, which has to be positive thus giving rise to a initial state 
without singularity, and (3 is the integration constant that contains the information 
about the initial velocity of expansion of the universe; then by recalling that r = 



a c we can use equation (3.23) to obtain the evolution of the shape of the universe 
as in the following 



a 
c 



D 



sx 



(3.26) 



in which we see that for t that tends to infinity the two factors tend to become 
constantly proportional, giving isotropization: by setting then D = 1 and choosing 
for simplicity 8X = —3mC\/ /3 2 — b the evolution of the single factors can be given 
as in the following 



3mC 



t + P-^^-b) 3 (b + 2pt + t 2 ) 



3mC 



t + P-y/0*-b 



(3.27a) 



(3.27b) 



with constraint given by 3m 2 (P 2 —b) = 1 for consistency with the fact that the time- 
time gravitational field equation is not actually a field equation but a constraint. 
Finally the spinor evolves as 



^3m(b+2/3t+t 2 ) 



( 



\ 







\ 



K 

C 



+ !)< 



-mt+\ 



■ In 



0) 



-mtH 



■ In 



!)) 







(3.28) 



/ 



with the spinor bilinears diluting as the inverse of the volume, correspondingly to 
the fact that they represent densities. By reading this solution we have that: at 
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the first instant of the evolution the universe has no singularities as it occupies 
a finite volume and it is shaped as a disc whose shorter axis is along the third 
coordinate described by the c factor, and it is filled with spinor fields with torsional 
interactions responsible for the anisotropy of the universe; as the time goes by the 
spinor field dilutes down, allowing isotropization of the universe without relenting 

2 

its expansion; eventually for very late epochs the scale factors have the ts scaling 
we know for the standard FLRW cosmic evolution. 

Now that we have the exact solutions of the f(R) = R case, we can employ them 
as the basis upon which to build the correction in 5a for the more general solutions 
of the f(R) = R + 5aR 2 . 

The form of the Dirac field will be unchanged: its form is going to give rise to a 
correspondingly simple form of the ip function given by 

<p = l- 5a— (3.29) 
r 

with potential 

V(<p) = Sa^r- (3.30) 

and the gravitational evolution equation for the volume of the universe is 

3mC\ „ mC ( t 2 t 3mC 
t 2 t At J 



( Z2 ~ ^ ~ — ) = ( 3 - 31 ) 



up to higher-order powers of 5a. As before, it is possible to integrate it to give the 
evolution of the volume of the universe as 



_ 3mC 



(b+2pt+t 2 ) + 5a^f W3 (P+t) In ( ^|±g±l )) (3.32) 

where q and £ are two new integration constants for which the initial volume of the 
universe is now different, but of course there is still no singularity: again these two 
new integration constants are linked by a relationship coming from the fact that 
the time-time gravitational field equation is a constraint. We remark that at later 
cosmological times in the evolution of the universe we still have isotropization, but 
more in general we have that this solution tends to be approximated to the same 
we would have had in the Einstein theory in what could be called Einsteinization. 



4 From the Jordan to the Einstein frame 

In this section we shall show that, by passing from the Jordan (gij) to the Einstein 
frame (gij) through a conformal transformation of the kind = ipgij (if <p > 0) 
or gij = —ipgij (if ip < 0), it is possible to quadratically integrate the Einstein-like 



equations (3.9) also in the general case f(R) / R. 



To start with, supposing for simplicity p > 0, we note that the relations between 
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the components of the tetrad fields associated with the metric tensors gij and gij 
are expressed as 

ef = v^ef , ^ = -jLeJ, (4.1) 



while the relation between the Levi-Civita connections induced by the two different 
frames is given by 



r 



Tij + 2<p dxi " l 2<p dxP 



(4.2) 



Moreover, we have also the identities 



(4.3) 



between the matrices V and f associated with the tetrads e and e respectively. 
Our aim is now to express the Dirac and Einstein-like equations in terms of the 
conformally transformed metric gij and its tetrad field e?. The covariant deriva- 
tives of the spinor fields are defined as in (3.6 ); in order to correctly evaluate them 
in the Einstein frame, we need to express the coefficients using the metric g~ij 
and its related quantities. With respect to this point, taking equations (2.3), (4.1 ), 
(4.2) and (4.3) as well as the identity 



(4.4) 



into account, we get the representation 



a 



i 



2<p 



2<p 



9jh [ r ik 3 - etae£ - ^d k <pSf + ^-d p <pgrig ik ) T h T k 



(4.5) 



Now, supposing that the metric gij is of the Bianchi type-I form (the metric g^ is 



then of the Bianchi type-I only up to a time reparametrization dt 



--dt) 



ds 2 = dt 2 - d 2 {t) dx 2 - b 2 (t)dy 2 - c 2 (t) dz l 



(4.6) 



the quantities (4.5) assume the explicit expression 



^o = 0, 


tlx = 


- - 


4^° : 








n 2 


= \h 2 i Q - 




= n 2 - 






n 3 


= ^C7 3 7° - 


- 




- ^ 



(4.7) 



Inserting equations (4.3) and (4.7) into equations (2.16), we obtain the following 
representation of the Dirac equations 



3<£> 



im 



3/ 



[(^)7° + i^VVVhV] i> = (4.8a) 
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j T 7 3(3 7 im - n 
* 2t Ap 



3i 



(4.8b) 

where we have defined f := a6c. Proceeding as in the Jordan frame, from equations 
(4.8) we derive the differential equations 

d 



dt 



[renvoi + —j{iipi b ip){^i°i>) = o 

3 



ar a/9? 



-r(^)(V'7 5 7V) = 



ife^ 7 5 7V)] + ^^(#7 5 V>) = o 



(4.9a) 

(4.9b) 
(4.9c) 



where Q := p 2 . Finally, equations (4.9) yield the relation (analogous of equation 

C 2 



(pi) 



,2 2 



(4.10) 



C being a constant. 

Now, we deal with the Einstein-like equations (2.15). As it has been shown in 
[HI El H21 [23 E], equations (2.15) can be expressed in terms of the metric Qij, 
assuming the form 



R, 



,, - \tigij = ^tf> + -^Sg - ^V&)gij + 7^3 (^757^)^757^)^ (4.11) 



where -Rjj and i? denote the Ricci tensor and the Ricci scalar curvature induced 
by the metric g~ij. Of course, also the tensors and E^ have to be represented 
making use of the metric §ij and the associated tetrad ef. In this regard, we have 



^fj '■= -(p + p)UiUj - -gij 
J p p 



(4.12) 



with UiUjg l i = 1. Concerning the tensor E^, from a direct calculation we obtain 
the identities 



J 00 



2c£ 



mipip 



16p 



(^757^^) (^757^^) 



E^ A = A = 1,2, 3 (4.13) 



while, due to equations (4.11), the vanishing of the non diagonal part of E^ gives 
rise to constraints analogous to those existing in the Jordan frame, namely 



J 12 



sg = o 

E? 3 = 



ab — 


ba 


= 


U 







cb- 


U 


= 


u 


'ipj 5 'j 1 ip = 





dc — 


ca 


= 


u 








(4.14a) 

(4.14b) 
(4.14c) 
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As in the Jordan frame, it is also seen that tfi A = (A = 1,2,3) are trivial 
identities, yielding no restrictions. 

Collecting all the obtained results, we conclude that in the Einstein frame (4.6) 

ssed as 

(V>757 T VO(V r )'57T'0) - -Kv(ip) (4.15a) 



the Einstein-like equations (2.15) are expressed as 

3 



a b be ac 

-7 + 7--+-- 

ab be ac 



P • 1 7/ 



2 2ip 2 rr 64(p 3 



b c be 

b c be 

a c ac 

- + - + -- 

a c ac 

a b ab 

a b ab 



P_ 



p_ 



1 
1 

1 

if 3 



V(<p) + ^3(^757 T V')(V , 757r'0) 

3 

V ( < P) + ^7-3(V'757 T V')(^757tV') 



64<p d 



V{<p) + 



64^- 



-(V , 757 t V ; )(V'757t^) 



(4.15b) 
(4.15c) 

(4.15d) 



together with the constraints (4.14). Again, the only significant case is when 
two components among a(t), b(t) and c{t) are equal. For instance, we suppose 
a{t) = b{t) and thus ^^^ip = ip^^tp = 0. We can now proceed as in the Jordan 
frame (see section 3), so obtaining from equations (4.15) the relations 

1 



D x exp X 



-dt 



(4.16) 
(4.17) 

The last equations we have to express in terms of the metric g~ij are those concerning 
the conservation law for the perfect fluid 



2-- :!V3^ + lm#4%) 



with 



(4.18) 
(4.19) 



Taking equations (|4.2|) and (|4.6|) into account, it is easily seen that equations (4.18) 

3 (p 



and (4.19) reduce to 



p - + d + xf- 

P r 



2(p 



;i + a) = o 



(4.20) 



Setting ® := tp a, one has — |^ = ^ and we can integrate equation ( 4.20| ) as 



p = po(t@) ( 1+a ) p = constant 



(4.21) 



At this point, we notice that if we could express the quantities p, t/jip e ip (or, 

3 

equivalently, G = as functions of f and insert the so obtained results into 

(4.22) 



equation (4.17), we would get a final equation for r of the kind 

r = /(f) 
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which is quadratically integrable. On the other hand, there exists a special though 
not too restrictive case for which the above mentioned conditions are satisfied. 
Indeed, setting ii/jj^ip = ^j 5 ^^ = in (4.10), we obtain 



C_ 
fQ 



equation (4.23) together with (4.21) and G 



(4.23) 



f'(R(ipip, p)) constitute a set 



of three relations involving the four variables r, ipip, p e 0. In principle, it is then 
possible to express the last three variables as suitable functions of f alone, thus 



obtaining from (4.17) a final quadratically integrable equation of the kind (4.22). 



As a last remark, we note that equation (4.15a) has the nature of a constraint 



on the initial data. As made in the Jordan frame, we have again to check that 
such a constraint is satisfied for all time after the initial instant. In connection 



with this, in Section 5 it is proved that the conservation laws (2.2) ensure that the 



four-divergence (with respect to the Levi-Civita connection associated with the 
conformal metric of the effective stress-energy tensor 



1 



1 



yD , 



1 3 

-^V^gij + (^757^X^757^)^ (4-24) 



vanishes. Repeating the arguments developed in section 3, we can apply again 
the Bruhat's result [16] and thus verify that the Hamiltonian constraint (4.15a) is 
preserved in time. 

We now illustrate two explicit examples where the above procedure actually works: 



Example 1: f(R) = R 4 model. To begin with, let us consider the model 
f(R) = R A coupled with a Dirac field (without perfect fluid for simplicity). In 
such a case, the corresponding potential (|2.8|) assumes the form 



V(<p) 



with the scalar field (2.7) given by 



1\ 5 3 7 

— (p3 



r' ~ I ( ^tp) 4 



(4.25) 



(4.26) 



Through equation (4.23) and the definition of G = tp we get the relation 

9 



Q = A 



C 



with 



.4 



, 3 

m\ 4 
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(4.27) 



(4.28) 
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from which, using again (4.23), we have 



3m 



32C 8 



and 



6 

tp 3 



9 iy 

4 V 4 / V4AC 8 



(4.29) 



(4.30) 



Inserting equations (4.29) and (4.30) into (4.17) and choosing the constant C 

(4.31) 



sufficiently small, we obtain the final differential equation 

f = Bf- 3 

with B > 0. The latter can be easily integrated as 



r = \ D{t + E) 2 + 



B 
D 



(4.32) 



where the constants D > and E have to be determined consistently with equation 



(4.15a). From equations (4.31) and (4.32) we have also r 2 > 5. We see that, at 



least in the Einstein frame, the model f(R) = R 4 coupled with just one Dirac Field 
gives rise to a universe with no initial singularity but undergoing a slow expansion 
without isotropization. 

Example 2: f(R) = R + aR 2 model. The second example we consider is given 
by the model f(R) = R + ctR 2 with a < 0, still coupled with a Dirac field alone. 
The associated potential (2.8) is now of the form 



V(<p) 



1 

8a 



while the scalar field (2.7) is 



(p = 1 — armjjijj 



(4.33) 



(4.34) 



From equation (4.23) and the definition of we get the algebraic equation of third 

(f 9 - amCf - f 3 9 = (4.35) 



order (for the unknown Gr) 



linking G to the volume-scale r. We recall that we are working under the assump- 
tion ip > 0. We are then forced to require Q > 0. In order to satisfy this condition, 



we have to impose that equation (4.35) has three real solutions (indeed, if equa- 



tion (4.35) had only one real solution, the latter would be negative as it is easily 



f > 



3^/3 



seen solving (4.35) by a graph). The imposed requirement yields the condition 



a\mC which represents a lower bound for the volume-scale r. Therefore, 
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again we can avoid initial singularity. Anyway, a solution remaining positive in 
time is 



2 r- 1 
Or = amC H — v3rcos I -< 

3 V 3 



where 



arccos 



3 y/3 \a\ mC 



(4.36) 



(4.37) 



At this point, choosing parameters such that a = — \ and mC = 1 for simplicity, 



end up with the final differential equation 



and inserting equations (|4.35|) and (|4.36|) into (|4.17|), after some calculation we 

'3\/3 N 



V3 



1 



3y/3 



cos - arccos , 
4 I 3 V 4f 



3^3 / 1 

H cos - arccos , 

16 I I 3 \ 4r 



(4.38) 



Equation (4.38) is quadratically integrable, but explicit calculations are not so 



easy. We can then proceed to numerically plotting it. For instance, if initial 
data f(0) = and r(0) = 1 are assigned, the corresponding behavior of the 
volume-scale f(t) during time is given by the graph in Figure 1 showing an (at 



3.x 10 s - 



2.x 10?- 



1.x 




2O0OO 



40000 60000 
t 



soooo 



100000 



Figure 1: volume-scale r vs time t 

least initially) accelerated expansion of the volume-scale r. Moreover, by plotting 
the function ?, it is seen that the volume-scale f in infinite of order n > 1. Due 



to equation (4.16), this fact implies that the universe isotropizes at large t. Also, 
we notice that at large t one has ~ 1 and the Einstein and Jordan frames tend 
to coincide. 
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5 Conservation laws 



We clarify the relationship between the conservation laws (2.2) and the Einstein- 
like equations (2.15 ). More in particular, we show that the validity the conservation 
laws (2.2) implies that the four-divergence of the Einstein- like equations (2.15), 
with respect to the Levi-Civita connection, vanishes and vice-versa. Also, we 
prove that the same result holds in the Einstein-frame, namely the four-divergence 
of the conformally transformed Einstein-like equations (4.11), with respect to the 
Levi-Civita connection associated with the conformal metric <?y = (pgij, is zero. 
To start with, for convenience of the reader, we rewrite equations (2.2) 



Vi&i + Ti&i - E ih T^ ih - l -S hiq R h ^ = 



v gijh + T oijh + S ij _ E ji 







(5.1a) 
(5.1b) 



where S JJ denotes the sum of the stress-energy tensors of the Dirac field and the 
perfect fluid. The curvature tensor can be decomposed in the form 



R 



h _ r>h , Y7 . 7v- h 

iqj ^ iqj ' v 3 qi 

h : 



V K- h +K P K h -K p K h 



where the contorsion tensor is expressed as the sum 11] 



ir.h i q, h 



with 



S 



K 



ij 



_ rr h 

f^ + f p g ph giJ 

1 dip 

2ip dxi 



(5.2) 

(5.3) 

(5.4a) 
(5.4b) 
(5.4c) 



Taking equations (5.2), (5.3) and (5.4) into account, a direct calculation shows 
that 



— -Si,- R hiq i 

2 ^hiq 1 1 



s h ip s\ p f h - \pW (s qth s^ h ) 



pV t [S h ^S j qh ) + <pVt (s^) s j 



qh 



(5.5) 



Moreover, making use of equations (|2.12|) and (|2.13|), it is easily seen that 

ViE ij + •/;>:'•' - r., h i jii ' 



V iS fe') + v iS M _ K^TP* - T j ih %W (5. 6 ) 

In |llj it is proved that = £W) — (pS hip S J ph and also that the conservation 

laws (5.1b I amount to the identities — £™ + V/ l S , - ?ife . From this and equations 
(5.3) and (5.4) one gets 



Vi£ W) = ViS (ii) - V 4 ( ipS hip S 



ph 



(5.7a) 
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K j^m = - v f h v q s h ^ + <ps ih j v q s hi(1 

- T^E^l = ipfiVqS™ - 2tpS jih V g S hiq 



(5.7b) 
(5.7c) 
(5.7d) 



where for simplicity we have defined ipi := J^. Inserting equations (5.5), (5.6) 



and (5.7) into (5.1a), the latter reduce to 



V; 



E(ii) _ IpV' f = 



(5.8) 



At this point, we recall that equations (2.7) and (2.8) are equivalent to the relation 

6. 



E VUp) + 2V'Up) = 



(5.9) 



(see [6] for the proof). After that, taking the trace of equation (2.15) as well as 



the identity -\S hqp S hqp = (V ; 757 t V ; )(V'757tV') P] into account, we get 



64</j 2 

-<pR - I-VW 1 + 3Vi^ + 4 V&) + 2<pS hgp S h iP 
lip ip 



(5.10) 



Substituting the identity E = E — ip§h qp S hqp [llj and equation (5.10) in equation 
(5.9), we obtain 



f, 3 1 i 

2 ip z 



-Vnp* + \v{ip) - -V'iip) - s hqp s hqp = 



(5.11) 



We rewrite equation (2.15) in the form 

1 



ipR 



" 2 ^^' : > 



7I1 



(5.12) 



The covariant divergence of (5.12) yields 



3 3 

■-<Pi(pj + -p h ^p h 9ij ~ V{p)gij 



Us hqp s hqp 



(5.13) 



By definition, the Einstein and the Ricci tensors satisfy V J Gjj = and {V 3 (p)R 



( V J VjVi — VjV-'Vj I Then equation (5.13) reduces to 



(5.14) 
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Finally, making use of equation (5.10) it is easily seen that 
"1 

.<P 

and then equations (5.13) amount to the equations V-'S(jj) — jf^Vj ( Sh qp S hqp J = 0, 



--fifj + -<p h <p gij 



V[ip)g^\ -\piS hqp S hqp (5.15) 



clearly identical to (5.8). 
Now, we denote by 



Tn = -%■) - ^V(p)g. 



(5.16) 



the effective stress-energy tensor on the right hand side of equations (4.11). Indi- 



cating by Vj the covariant derivative associated with the conformal metric = 
tpgii, we show that the condition V-'Tjj = is equivalent to the conservation laws 

V J S (<i ) - \ipVi (s hqp S hqp ) = 0. Recalling equations gjg , we have 



ij ~ g v s-Lii — g 

<p p 



» ' sTij 



2v? s <9x s 4 



1 



2p V <9x? S <9x s 



(5.17) 



We have separately 



V'V S T^ = Vv s - ^V((p)gi3 - \s hqp S hqp g t ^ 

V{<p) 



l«7- (i + V 



1 „„- 1 / dp r0 , dip q 



i s 3 - 



<p> lip \dx 



2<p 
dp 



(5.18) 



dx 



s i 



dx u 



g uq g S i T. 



- 13 



-9 



s 3 . 



p 2<p \dx 



dx s 



dip 
dx u 



g uq g si 



ip- 



2p 



1 g*i ( 



, dip - ~ (u 



~Ap^ gS:> {^x~i^ hqp ^ hqP9sj + ^^hqpS^gij - -7^-S hqp S hqp g si 



+ 



(5.19) 



1 dip - 2 dip 1 dyw 

2^3 dx* w c^ 2 <9x* h9P 
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and 



</? 2y? V^ J ' ' ' n "* ' M -" 



/l - 1 

x o 



1 



3^ Sh qp S hqp gis + Sh qP S hqp gij S hqp S hgp g sj 



X 



+ (5.20) 



2<f?dx* ah9pL 

Collecting equations ( |5.18| ), ( |5.19| | and ( |5.20[ ) and recalling that E = t-ipS hqp S hqp , 
we have then 



2 (j9 



1 



<(ij) 2 ^i^hqpS qP 



(5.21) 



because the identity — |E + ^V(¥0 — V'Cy) = holds identically, being equivalent 
to the definition 99 = f{F(T)) [6]. 



6 Conclusion 



In the present paper we started from gravitational theories of the f(R) type in 
which the Ricci scalar R contains torsion as well as the metric field, considering 
them in the case in which the spacetime is filled with Dirac matter fields; we have 
employed these models to study anisotropic cosmological models BI: with respect 
to Einstein gravity, additional gravitational effects arise as a consequence of the 
non-linearity of the f(R) function and the natural torsion-spin coupling, known to 
induce centrifugal barriers [22]. These non- linear and repulsive centrifugal effects 
may enforce one another, especially in the case anisotropies are considered; the 
relationship between torsion and anisotropies of the spacetime has been widely 
investigated in the literature, and recent accounts are for example those in |23U24j ). 
An important issue that must be highlighted is the fact that, despite the anisotropic 
background, the Einstein tensor is diagonal, while, because of the intrinsic features 
of the spinor field, the energy tensor is not diagonal: in this circumstance the non- 



diagonal part of the gravitational field equations results into the constraints (3.10) 
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characterizing the structure of the spacetime or the helicity of the spinor field or 
both; in our understanding, the only physically meaningful situation is the one in 
which two axes are equal and one spatial component of the axial vector torsion 
does not vanish, giving rise to a universe that is spatially shaped as an ellipsoid of 
rotation revolving about the only axis along which the spin density is not equal to 
zero. It is also worth noticing that for the specific form of the Dirac field we have 
chosen, although the Dirac field is massive, still we have the conservation of the 
axial current; on the other hand it is widely known that quantum corrections might 
produce additional terms for which the axial current is not conserved anymore, if 
the Dirac field is charged or if additional terms in the dynamics are introduced, 
and then it would be interesting to see what are the effects that constraints such 
as (3.10) have on these partially conserved axial currents. However because in the 
present paper we do not deal with quantum corrections, nor with charged fields or 
extensions in the electrodynamics, addressing this issue would bring us far from 
the main aim of this work, and we refer the reader to [18] and references therein 
for a discussion on this subject. 

The field equations have been worked out in both the Jordan and the Einstein 
frames, showing that the procedure to get solutions proposed in O [To] applies 
also here, independently of the given function f(R). In particular, in the Einstein 
frame the Einstein-like equations result in general to be quadratically integrable. 
As it may have been expected, these models are intrinsically quite difficult to study 
in general; we have focused on some specific examples for which the field equations 
are considerably simplified: we have studied the model f(R) = -R 4 , where exact 
analytic solutions are found, although this model does not reproduce the physical 
content we expect; we have further studied the model f(R) = R + aR 2 , for which 
exact analytic solutions are very difficult to obtain; therefore, in this case, we 
have first worked in the Jordan frame, finding solutions that are analytic but 
approximated to the first order of a, and then we have moved in the Einstein 
frame, employing numerical methods. In the case f(R) = R 4 , we have found that 
the universe does not display isotropization, as it should, and there is no correct 
Einstein limit; in the f(R) = R+aR? case, we have found that the results obtained 
in the two frames coincide, and they yield Einsteinization and isotropization, that 
is as the cosmic time goes by both the non-linearities of the f(R) function and the 
anisotropies of the background tend to vanish. In both models however, no initial 
singularity for the volume of the universe has been found at all. 
Finally, we have proved that the proposed method for solving the field equations is 
consistent with the assignment of initial data, clarifying the role of the previously 
stated general conservation laws [TT] in preserving the Hamiltonian constraint. 
To summarize, we have found here three types of results: first of all, from the point 
of view of the behaviour of the universe, our results may help to solve one of the 
major conceptual problems of f(R) theories related to the arbitrariness of the f(R) 
function; indeed, as we have shown, a given f(R) may either give rise to physics 
with no good behaviour or yield physics with good behaviour; thus although there 
is no theoretical reason to prefer a particular form for the f(R) function instead of 
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another, nevertheless there are phenomenological set-ups in which selection rules 
discriminating physical f(R) function from infeasible ones may be established. 
Secondly, universes that have finite initial volume are possible, and therefore the 
gravitational singularity problem receives new fuel from the present discussion. 
Finally, the initial values problem is shown to be well formulated. 
The general outlook emerging from our study encourages to pursue the study of 
torsional f{R) models, since their field equations and conservation laws are well 
defined and consistent with the assignment of initial data, while the presence of 
torsion can remove the singularity problem, and there are hopes that a possible 
specific f(R) may be selected which will be able to describe what is missing in the 
evolution of the universe as we know it. 
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